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1 Introduction 

Recently, nucleus-nucleus collisions at high energy have been studied very inten- 
sively because of the possibility of the formation of a quark-gluon plasma state. Also, 
the rapidity distributions of the charged particles and the ratios of the numbers of 
anti- hadrons to hadrons at SPS[1] and RHIC[2]. 

The nucleon-nucleus(A r — A) and nucleus-nucleus(^ — A) collisions have been stud- 
ied from the multiple-scattering view[3,4,5] . In particular, the multi-chain model(MCM) 
succeeded phenomenologically to reproduce the old data on the multiplicity and the 
inclusive spectra of the leading nucleon and secondary particles etc. in N — A and 
A — A collisions[6,7,8,9]. In MCM[10], it is assumed that the nucleon-nucleon(A f — N) 
interaction is considered as a exchange of one multi peripheral chain from which sec- 
ondary hadrons are emitted. However, it was difficult to perform analytical calculation 
for A — A collisions because of the complexity of the general formula in resolving the 
full combination of chain configuration. Thus, dynamical simulation models for A — A 
collisions have been developed. 

In this paper, on the basis of MCM with successive collision[ll], the full combi- 
nations of chain configuration are resolved in terms of the vector-operator method in 
the moment space by consulting the Mellin transformation[12]. The analytic forms for 
single-particle distributions in A — A collisions are derived. 
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In this paper, the dynamics is assumed to scale with energy. A nucleus is treated 
as a set of mutually independent nucleons. The cascading of the produced hadrons is 
neglected owing to the long formation length in nucleus. Also, the longitudinal motion 
is treated. The transverse momentum distribution is assumed to be independent of 
incident energy and longitudinal momentum. 

In Section 2, the MCM with successive collision for N — A collisions is summarized 
and the vector-operator formalism in the moment space is given. In Section 3, single- 
particle distributions are investigated and their analytic forms are derived. In Section 
4, conclusion and discussion are given. 



2 MCM with successive collision for N — A col- 
lisions 

2.1 Structure of the inelastic N — N interaction 

We start with the inelastic N — N interaction. We assume that the inelastic N — N 
interaction as the exchange of one chain with a mesonic cluster from which secondary 
hadrons are emitted for simplicity. Namely, we consider the reaction 

N + N^N + N + M^N + N + hadrons 

where M denotes the mesonic cluster. We pay attention to the single-particle distri- 
bution of the inclusive process N + N — ¥ C + A (anything) given by 

NN da(NN -+CX) 
Pc (x) = x — (x > 0) (1) 

where x denotes the longitudinal momentum fraction(Feynman variable). We express 
the projectile(target) fragmentation regions asa;>0(x<0). 

By the Mellin transformation of Eq.(l), we define the moment Pc N (J) and the 
operator p NN (J) in the moment space as 

Pc N {J)= f dxx J - 2 p^ N {x)^<C\p NN (J)\N> (2) 
Jo 

where |7V > denotes the nucleon state in the projectile nucleus A. Here, we employ a 
vector-operator notation[12]. We introduce the operator J(J) for the inelastic inter- 
action of N — Acollisions as 

J(J) = ~iwP NN (J) (3) 

inel 

where <t^^J is the inelastic cross section of N — N collisions. 

We define the matrix elements of the leading nucleon N which is not newly emitted 
from the mesonic cluster and the mesonic cluster M as 

< N\J(J)\N >= -Lf < N\p NN (J)\N > = F(J), (4) 

< M\j(J)\N >=^< M\p NN (J)\N K(J) (5) 

inel 
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where < N\ and < M\ denote the leading nucleon state and the mesonic cluster state, 
respectively. By the inverse Mellin transformation, the fraction functions F(x) and 
K(x) are given by 



2 « L- 



F( x ) = 7^7 / dJx'-'FiJ), (6) 



C+lOO 

1-J j 



K{x) = ^- I dJx'-'K(J) (7) 
where f ^F(x) — 1 and K(x) = j^F(l — x) from the cascade mechanism. 

2.2 MCM with successive collision and recurrence equa- 
tions 

In MCM with successive collision for N — A collisions, the basic assumption is that 
only the projectile nucleon can interact successively with the nucleons inside the target 
nucleus A, while nucleons inside nucleus A experience only one inelastic collision with 
projectile nucleon. This model can reproduce the single-particle distributions of the 
leading nucleon and secondary particles. We show a sketch of MCM in Fig.l. 

First, according to Ref.[ll], we summarize the MCM with successive collision for 
N — A collisions at the impact parameter b. The inclusive distributions of leading 
nucleon and mesonic cluster after n-times successive collisions of the projectile nucleon 
with nucleon inside the target nucleus A in the projectile fragmentation regions(a; > 0) 
are given by Q(x, n; b) and M(x, n; b), respectively. They satisfy the integral equations 



C^fc-) 



Q(x,n;b) = (l-X(b)Q(x,n-l;b)+X(b) -»F(-)Q{y,n - 1; 6), (8) 



M(x,n;b) = M(x,n-l;b) + X(b) [ ^-K(-)Q(y, n - 1; b) (9) 

Jx y y 

where the functions F(x/y) and K(x/y) are given by Eqs.(6) and (7). Also, X(b) 
is the inelastic interaction probability defined as X(b) = crf^TAib) where Ta(o) is 
the nuclear thickness function of nucleus A normalized to unity. It is noted that the 
collision number n is limited to < n < A where A denotes the mass number of the 
nucleus A. 

By the iteration method, we obtain the following solutions of Eqs.(8) and (9) in 
the final collision^ = ^4): 



Q(x, A; b) = P (A; b)S(l -x) + J2 Pm{A; b)F (m) (x), (10) 

A m— 1 „ i 

M(x,A;b) = J2 P m(A;b)J2 ^-K&F (l) (y) (11) 

m=l 1=0 Jx V V 
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Fig.l MCM with one kind of chain for N — A collisions. Wavy lines represent the 
inelastic N — N interaction. 

where F^\x) = £ ^F{'±)F^-^ (y) and F^(y) = 8(1 ~ y). Also, P m (A;b) = 

m^j A(b)) A ~ m A(b) m which is the Glauber probability. The inelastic cross section 

of N — A collisions is given by of^i = J d^X)m=i Prn{A; b) and the averaged collision 
number n — Aa^i/viL^i- The moments of Eqs.(lO) and (11) are given by 

A 

Q(J, A; b) = P (A; b) + J2 P ™( A '' b)F(J) m , (12) 

m — 1 

A m— 1 

M(J,A;b) = ^P m (A;b)^^(J)F(J) i . (13) 

m=l 1=0 

Next, we reformulate the above MCM with successive collision in order to extend 
the model to A — A collisions. 
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(i) The operator Qa{ J, n; b) is introduced and the moments Q(J, n; b) and M(J, n; 6) 
are defined as follows: 

Q(J,n;b)=<N\Q A {J,n;b)\N>, (14) 



M(J,n;b) =< M\Q A (J,n;b)\N > . (15) 

We assume the following recurrence equations for Q(J,n;b) on the analogy with 
Eq.(8) in the moment space: 

Qa(J, n; b) = [(1 - \{b)) + \{b)J(J)]Q A (J, n - 1; b) (16) 

where Q(J,n; b) = 1. From Eq.(16), we get the solution 

Qa(J, n; b) = [(1 - \(b)) + \(b)j(J)] n . (17) 

Then, we put n — A and use the binomial expansion. Eq.(17) reduces to 

A 

Q A (J,A;b) = P (A;b) + J2 P ^;b)j(Jr. (18) 

m — 1 

where m agrees with the number of chain in MCM. 

(ii) We define the matrix elements for J(J) m in Eq.(18) as 

N(J,m) =< N\j(J) m \N >, (19) 

M(J,m) =< M\j(J) m \N > (20) 

where N(J, 0) = 1 and M(J,0) = 0. Also, N(J, 1) = F(J) and M( J, 1) = K(J) from 
Eqs.(4) and (5). 

We assume the recurrence equations 

N(J,m) = F(J)N(J,m-V), (21) 

M( J, m) = M(J, m - 1) + K(J)N(J, m - 1). (22) 

Eq.(21) satisfies the baryon number conservation and Eq.(22) is derived on the analogy 
with Eq.(9) in the moment space. The solutions of Eqs.(21) and (22) are given by 

N(J,m) = F(jr, (23) 



M(J,m) = ^2K(J)F(J) 1 . (24) 

1=0 

Therefore, from Eqs.(18) and (23), Eq.(12) is derived and from Eqs.(18) and (24), 
Eq.(13) is derived. Thus, Eqs.(10) and (11) of MCM with successive collision for 
N — A collisions are derived from the three recurrence equations of Eqs.(16), (21) and 
(22). 
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3 Single-particle distributions in A — A collisions 



We investigate single-particle distributions in A — A collisions on the basis of 
Eqs.(16), (21) and (22) for TV — A collisions. Consider the collision of two nuclei 
with mass numbers A and B(A — B collisions) at the impact parameter b. We treat 
the inclusive processes of AB — ¥ NX and MX for the projectile fragmentation regions 
of nucleus A(x > 0). 

First, we introduce the operator Qab(J; b, A) in the moment space. Consider the 
moments Qj^ N (J; b, A), M A 4 B N (J;b, A) of the inclusive distributions Q A , B N {x;b, A) 
and M AB N (x; b, A). They are defined as 

Qif N (J; b, A) =< N\Q AB (J; b, A)\N >, (25) 

M$ N (J;b,A) =< M\Q AB (J;b,A)\N > (26) 

where the argument A denotes the treatment of particles observed in the fragmentation 
regions of nucleus A in final state. The \N > is the nucleon state inside the projectile 
nucleus A and the states < N\ and < M\ denote the leading nucleon and the mesonic 
cluster observed in the projectile fragmentation regions of nucleus A in the initial state, 
respectively. 

Second, we introduce the operator Qab{J, n; b) which characterizes the total struc- 
ture of A — A collisions. The recurrence equation is assumed as follows 

Q AB (J,n;b) = [(l-X{b))G(J) + \(b)J(J)]Q A B(J,n-l;b), (27) 

where n is the collision number(0 < n < AB) and Qab(J, 0; 6) = 1. Also, A(6) = 
Pitiel fdsf dtT A {s)T B {t)S 2 (b + s-t) where T A (s) and T B {t) are the normalized nu- 
clear thickness functions of nucleus A and nucleus B, respectively. We introduce the 
operator G(J) so as to express the passing-through probability (1 — A(6)) explicitly. 
By considering both operators G(J) and J (J), we may calculate full combinations of 
the chain configurations with ease. In fact, it corresponds to MCM with two kinds of 
chain as shown in Fig. 2. 

The solution of Eq.(27) is given by 

Q AB (J,n;b) = [(l-\(b))G(J) + \(b)J(J)} n . (28) 

The operator Q A b (J; b) is defined to the limit 

Q AB {J;b)= lim Q AB (J, n; b). (29) 

Taking notice of the ordering of G(J) and J(J), we expand Q A b(J;1>). Then, 
Eq.(29) reduces to 

Q AB (J;b) = (1 - \(b)) AB G(J) AB + (1 - A(&)) AS - 1 A(&)[G(J) AS - 1 J(J) 
+ G{J) AB - 2 J{J)G{J) + ■■■ + Ji^GiJ)^- 1 ] + ■■■ + \{b) AB J(J) AB . (30) 



G 



A 




Fig.2 



Fig. 2 MCM with two kinds of chain for A — B collisions. Wavy lines and dashed 
lines represent the inelastic interaction and the passing-through in N — N collision, 
respectively. 

If G(J) = 1, Eq.(30) is given by 

AB 

Q(J;b) = P (AB-b) + P m {AB-b)J{J) m (31) 

m—1 

where P m (AB; b) = (^^{1- \{b)) AB - m \{b) m that is the Glauber probability for 

A — B collisions and m agrees with the number of chain in MCM. The inelastic cross 
section of A — B collisions is given by = J rf&X)m=i Pm(AB : b) and the averaged 
collision number n = ABa^ l /crf^ l . 

Third, we discuss the relationship between Qab(J', b, A) and Qab{J', b). In MCM, 
the maximum number of collision of one nucleon inside the projectile nucleus A is 
restricted to the mass number of the target nucleus B as investigated in iV — A collisions 
and vice versa. Thus, we must impose the restriction to Qab{J\ b). When we observe 
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the particle in the projectile fragmentation regions, we introduce the rule of the sum 
decomposition in the projectile fragmentation regions of nucleus A for the operator 
product 0\Oi ■ ■ ■ Oab in Eq.(30) as follows: 

A 

: 6i0 2 - • -Oab :a= ^ nf =1 6 j+(i _ 1)B . (32) 

i=l 

iFrom Eq.(32), we get 

: Q{J;b) : A = (1 — \(b)) AB AG{J) B + (1 — A(b)) AS_1 A(b)[( J 4 — 1)G(J) B + G(J) B ~ 1 J(J) 

+ ■■■ + JiJW) 13 - 1 + (A- 1)G(J) B ] + ■■■ + \(b) AB AJ(J) B (33) 
Furthermore, we use the relation [G(J), J (J)] = 0. Eq.(33) reduces to 

AB B 

: Q(J; b) : A = A ^ Pm(AB; b) ^ H(l; AB, m, B)G{J) B - 1 J{J) 1 (34) 

m=0 1 = 

where H(l; AB,m, B) is the hypergeometric distribution defined as 



H(l;AB,m,B) 



AB -m (m 
B-l ) \ I 

AB 
B 



The operator Qas(J;&, A) is defined as 

Qab{J-XA) =: Q{J;b) : A . (35) 
Thus, from Eq.(34), we obtain 

AB B 

Q AB (J;b,A) = A^P m {AB;b)^H{l;AB,m,B)j{J) 1 (36) 

m=0 1=0 

where we put G(J) = 1. From Eqs.(25), (26) and (36), 

AB B 

Qirf N ( J; b, A) = AJ2 p m( A B;b)J2 H(l;AB, m,B)F(J)\ (37) 

m=0 1=0 



AB B l-l 

M^%(J; b,A) = A^2 Pm(AB; b) H ih AB, m, B) ^ K(J)F{J) r . (38) 

m=l 1=1 r=0 

By means of the inverse Mellin transformation of Eqs.(37) and (38), the single-particle 
distributions of AB — > NX and MX for the projectile fragmentation regions(a; > 0) 
are given by 

/AB B 
dbA^P m {AB : b)^H{l;AB,m,B)F { - l \x), (39) 

m=l (=1 
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/AB B l-l „! 

dbA^P m {AB:b)^H{l;AB,m,B)^2 I ^K{-)F (r \x). (40) 

m=l (=1 r=l"' a; ^ ^ 

We obtain the similar results for the target fragmentation regions(:r < 0). ft is easy to 
show that Eqs.(37) and (38) reduce to Eqs.(12) and (13) for x > in N — A collisions. 
Also, we obtain 

p% A (x) = A*ZeiF(x), (41) 

Pm A (x) = AaZeiK(x) (42) 
for the target fragmentation regions of nucleus A(x < 0). 

4 Conclusion and Discussion 

We have investigated the single-particle distribution in A — A collisions on the basis 
of MCM with successive collision. By assumption of three recurrence equations with 
the rule of the sum decomposition, we derive the analytic form of the single-particle 
distribution of the inclusive process A + B — ► C + X which is expressed in the Glauber 
probability with the hypergeometric distribution. 

This result may be interpreted from the probabilistic view as follows: The total 
number of collisions in A — B collisions is AB. The maximum number of collisions 
one nucleon inside projectile nucleon A is B. The number of the inelastic N — N 
interaction( J(J) ) and the number of the passing-through (G(J)) correspond to m 
and AB — m, respectively. If the number of the inelastic N — N interaction actually 
occurring in the event is I, the probability of occurring exactly I can be calculated by 
the hypergeometric distribution^^; AB, m, B)) . 

Recently, the polarization phenomena in A — A collisions are payed attention to 
and the global polarization[13] is investigated. This theme is one of the interesting 
problems in future. 



References 

[1] W. Bell et al., Z. Phys. C27, 191 (1985). 

[2] G. Bearden et al.,Phys. Rev. Lett. 88, 202301(2002); ibid. 90, 102301 (2003). 

[3] R. J. Glauber, Lectures on Theor. Phys. (Springer, Berlin) 1, 315 (1959). 

[4] R. J. Glauber and G. Mattehiae, Nucl. Phys. B21, 135 (1970). 

[5] A. Cappella and A. Krzwicki, Phys. Lett. 67, 84 (1977). 

[6] K. Kinoshita, A. Minaka and H. Sumiyoshi, Prog. Theor. Phys. 61, 165(1979). 

[7] S. Date, M. Gyulassy and H. Sumiyoshi, Prog. Theor. Phys. 63, 928 (1980). 

[8] K. Kinoshita, A. Minaka and H. Sumiyoshi, Z. Phys. C8, 205 (1981). 

[9] H. Sumiyoshi, Phys. Lett. B113, 241 (1983). 

[10] T. Kanki, K. Kinoshita, H. Sumiyoshi and F. Takagi, Prog. Theor. Phys. 
Supplement 97B, 134 (1989). 



9 



[11] H. Noda and T. Tashiro, Z. Phys. C45, 85 (1989). 

[12] K. Kinoshita, H. Noda, T. Tashiro and M. Mizouchi, Z. Phys. C4, 103 (1980). 
[13] Z. -T. Liang and X. -N. Wang, Phys. Rev. Lett. 91, 102301 (2005). 



10 



